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Abstract. We consider a system of independent branching random walks on 
R which start off a Poisson point process with intensity of the form e\(du) = 
e~^ u du, where A e 1 is chosen in such a way that the overall intensity of 
particles is preserved. Denote by \ the cluster distribution and let tp be the log- 
Laplace transform of the intensity of \. If > 0, we show that the system 
is persistent (stable) meaning that the point process formed by the particles 
in the n-th generation converges as n — > oo to a non-trivial point process Tl* x 
with intensity e\. If Ai//(A) < 0, then the branching population suffers local 
extinction meaning that the limiting point process is empty. We characterize 
(generally, non-stationary) point processes on R which are cluster-invariant 
with respect to the cluster distribution \ as mixtures of the point processes 
II*, A over c> and A G K st , where K st = {A e R : <p(X) = 0, Xip'(X) > 0}. 



1. Introduction 

1.1. Persistence criterion. We consider a population of particles on R d whose 
stochastic evolution (branching dynamics) is governed by the following rules. The 
initial positions of the particles form a Poisson point process on M. d denoted by tto. 
The intensity measure of tto, denoted by v, is assumed to be finite on bounded sets. 
After one unit of time each particle is replaced, independently of all other particles, 
by a random cluster of offsprings whose displacements with respect to the parent 
particle are governed by a point process \ (which is allowed to be empty) . Note that 
the distribution of \ does not depend on the position of the parent particle. The 
daughter particles form a point process denoted by m. Then, any of the daughter 
particles is independently replaced by a random cluster of granddaughter particles 
forming a point process 7T2, and so on. We denote the point process formed by the 
particles of the n-th generation by 7r n . We will be interested in the behavior of n n 
as n — > oo. 

This problem has been much studied especially if the cluster distribution \ is 
critical (i.e., the mean number of particles in the cluster x is 1) and v is the Lebesgue 
measure; see 14, 11, loL 21 | and [Lsl Chapters 11, 12]. In this particular case it has 
been shown, under a second moment condition on the intensity of the clusters, that 
in dimensions d — 1 and d — 2 the branching dynamics suffers local extinction. This 
means that the point process 7r„ converges as n — > oo to the empty point process. 
In dimension d > 3, the branching dynamics is persistent (stable) meaning that ir n 



converges to some non-trivial point process whose intensity is equal to the Lebesgue 
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measure. A general theory of convergence to equilibrium of branching populations 
on an arbitrary Polish space was developed in [lij]; see also po| . 

We will be interested what happens if the Poisson point process ttq formed by the 
initial positions of the particles is not homogeneous, i.e., the measure v is not the 
Lebesgue measure. We will need the following assumption. Let J be the intensity 
of the point process x, that is J is a measure on R d given by J(B) = E[\(B)] for 
every Borel set B C M. d . Assume that the log-Laplace transform of J is finite: 



ip(t) := log / e {t ' u) J(du) < oo, t £ 



(1) 



Here, (•, •) denotes the scalar product on R d . 

When studying persistence of a branching dynamics it is natural to require that 
the intensity of the point process 7r„ should remain constant in time. This leads to 
a convolution equation 

(2) J * v = v. 

The solutions of this equation have been described by Deny [9j]. For A £ R rf let e\ 
be a measure on R d whose density with respect to the Lebesgue measure is given 

by 

(3) e x (du) = e-( x > u >du. 

An easy computation shows that any measure of the form e\ , where A is such that 
f(X) = 0, is a solution to ([2]). Conversely, under a non-lattice assumption on J, 
Deny [9( showed that any measure v which solves ((2]) and is finite on bounded sets 
can be represented as a mixture of the form 



(4) 



v{-) = / e x (-MdX), 



K 



where fi is a finite measure on the set K := {A e M. d : (p(X) = 0}. 

In the next theorem we give a criterion for persistensy /local extinction of the 
branching dynamics starting off a Poisson process with intensity measure e\, AeX. 
Denote by V the gradient operator. 

Theorem 1. Consider a branching dynamics satisfying ^ which starts off a Pois- 
son point process with intensity e\, where A £ M. d is such that (p(X) = 0. 

(1) If (A, Vy?(A)) > 0, then the point process ir n converges weakly to some point 
process IT* having intensity e\ (the branching dynamics is persistent) . 

(2) If (A, Vy(A)) < 0, then the point process 7r„ converges weakly to the empty 
point process (the branching dynamics suffers local extinction) . 

1.2. The one-dimensional case. It is instructive to consider the special case 
d = 1. To exclude trivialities, we assume that J is not concentrated on [0, +oo) or 
(— co, 0]. This implies that limt^-too ip(t) — +oo. The function ip is strictly convex. 
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Assume first that the cluster distribution \ is subcritical, that is, the mean 
number of offsprings of any particle is less than 1. Note that this implies that the 
progeny of any individual particle dies out with probability 1; see |2j. Let us see 
what happens to the population as a whole. Since the subcriticality means that 
(p(0) < 0, the equation ip(X) = has two solutions Ai < and A2 > which satisfy 
V?'(Ai) < and ^'{M) > 0. By Theorem Q] both measures e\ 1 and e\ 2 give rise 
to a persistent branching dynamics. In both cases we have a system of branching 
random walks which become extinct individually but survive as a collective. 

Assume now that the cluster distribution \ is critical, that is the mean number 
of offsprings of any particle is 1. Equivalently, y(0) = 0. Again, the progeny of any 
individual particle dies out with probability 1; see [2j. If <p'(0) = L uJ (du) = 0, 
then the equation <p(\) ~ has only one solution A = 0. It is known, see [14]], that 
the corresponding branching dynamics, starting off the Lebesgue measure, suffers 
local extinction. Therefore, we assume that J R uJ(du) < which means that the 
particles have a drift to the left. (The case J R uJ(du) > is analogous.) Then, 
(f(0) = and tp'(O) < 0. Hence, the equation <p(\) — has two solutions Ai = and 
A2 > 0. By [l4(, the first solution Ai = gives rise to a branching dynamics which 
suffers local extinction. However, by Theorcm[IJ the second solution A2 corresponds 
to a persistent branching dynamics. Thus, a critical branching dynamics can survive 
provided that the particles have a non-zero drift and they have chosen their starting 
positions in a right way. 

Assume finally that the cluster distribution \ is supercritical, i.e., the mean 
number of offsprings of any particle is greater than 1. Equivalently, <ys(0) > 0. 
There are 3 possibilities. First, it may happen that the equation tp(\) = has no 
solutions, which is not interesting to us. Second, it may happen that the equation 
<p(\) = has two different solutions of the same sign. By Theorem [1] the solution 
which has larger distance to gives rise to a persistent branching dynamics, whereas 
the branching dynamics corresponding to the remaining solution becomes locally 
extinct. Thus, a supercritical branching dynamics may suffer local extinction, even 
though the progeny of any particle survives with positive probability and grows 
exponentially on survival. Finally, it may happen that the equation ip(X) = has 
only one solution Ai such that ip'(Xi) = 0. This boundary case is not covered by 
Theorem [TJ In Section 12.41 we will show that the branching dynamics dies out in 
the boundary case if the underlying branching process is the branching Brownian 
motion, and conjecture that this holds for general branching random walks. The 
boundary case has been an object of recent interest due to its connection to the 
extremal process of the branching Brownian motion; see [TH], @, HI, [H ■ 

A very simple special case arises if the cluster x consists a.s. of exactly 1 particle 
(no branching). The persistence of the system is evident, since 7r„ is a Poisson 
process with the same intensity e\ for every n. See [H K3l EH, [l9[ for various types 
of exponential intensity particle systems without branching. 

2. Properties of the competing branching dynamics 

In the one-dimensional case d = 1 the branching dynamics described in the 
introduction can be seen as a competition between branching particles. In the 
sequel, we will study the properties of this competition always assuming d = 1. 

2.1. Notation. Let M be the space of discrete point configurations (locally finite 
counting measures) on R. A point process is a random element of M. Given a 
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discrete point configuration p — Sv t S M, where the index set / is at most 

countable, we denote by max/? := maxig/ Vi the position of the right-most atom 
(the leader, the most successful particle) of p (provided that it exists). Also, for 
a£l we denote by T u p — X^p/ ^«+Vi ^ ne P om t configuration obtained by shifting 
every atom Vi of p by u. Let No = {0, 1, . . .} be the set of non-negative integers. 

We fix some point process \ (t ne distribution of clusters) with intensity J satis- 
fying ([I}. Note that the number of particles in x is a - s - finite. A branching random 
walk with cluster distribution x is a process which starts at time with a single 
particle located at the origin and in which every unit of time every particle is in- 
dependently replaced by a cluster of daughter particles whose displacements with 
respect to the parent particle are governed by the point process \- We denote by 
Xn the point process formed by the particles in the n-th generation of the branching 
random walk. By the Hammersley-Kingman-Biggins theorem 0, 0] , the maximum 
max x-n of a supercritical branching random walk Xn satisfies 

(5) lim — maxx„ = (3q a.s. on survival, 

where (3q is the largest zero of I(z) — sup teR (zi — ip(t)), the Legendre-Fenchel 
transform of ip. Note that / is allowed to take the value +oo. 

In the sequel, we denote by ttq = J>2iLi ^Ui a Poisson point process on R with 
intensity e\(du) — e~~ Xu du. With this notation, the positions of the particles in the 
n-th generation of a branching dynamics starting off ttq form a point process 

oo 

(6) ^ = E T ^ } )> 

1=1 

where {xn^ ■ n £ No}, i £ N, are copies of the branching random walk {%„ : n £ 
No}, all objects being independent. 

2.2. The leading particle. Let us look at the progeny of an individual particle in 
the branching dynamics {7r„ : n € No}. We claim that even in the persistent case, 
the progeny of any particle with probability 1 either dies out or drifts to +oo or 
— oo. Indeed, if the cluster distribution x is sub-critical or critical, then the progeny 
of any particle dies out a.s.; see 0. If the cluster distribution is supercritical, then 
the progeny survives with positive probability. However, assuming for concreteness 
that A > and recalling that <p(\) = 0, it is easy to see that the limit /3q in ([S]) is 
negative. Hence, the progeny of any individual particle either dies out or drifts to 
— oo in a uniform way. 

These considerations suggest that the leader of a persistent branching dynamics 
at time n, where n is large, is most likely to be an offspring of some particle which 
had a very unfavorable position at time (we assume A > here). The next 
proposition makes this precise. Let i n denote the (random) number of the particle 
which is the ancestor at time of the leader of the system at time n. That is, i n is 
such that max7r„ = Ui n + max^n • 

Proposition 1. Assume that ip(X) — 0, A > 0, <p'(X) > 0. Then, we have 
liron^oo \U ln = < a.s. 

In the supercritical or critical case, let /3o be the largest zero of the Legendre- 
Fenchel transform /. In the subcritical case, we set f3o — — oo. In the next proposi- 
tion we show that the maximum of a supercritical branching dynamics 7r n suffering 
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local extinction converges to — oo with the same asymptotic linear speed flo as the 
most successful offspring of a single particle. 

Proposition 2. Assume that ip(X) — 0. A > 0, <p'(X) < 0. Then, we have 
lirn n _ ! . 00 i max7r n = /3 a.s. 

Let us mention also the following proposition which will be obtained as a simple 
corollary of Theorem [T] 

Proposition 3. Let {%„ : n € No} be a branching random walk with cluster distri- 
bution x which satisfies (TT]). Let t > 0. 

(1) If(p'(t) > P , then c(t) := ]im n ^ 00 e -*'( t ) n E[e* miBt ^»] exists in (0,oo). 

(2) If iff it) < fa, then c(t) := Iim^oo e -^( f )™E[e tmaxx "] = 0. 

Remark 1. We conjecture that for (f'(t) — /3q the limit is 0. The proof in the case 
of the binary branching Brownian motion will be given below; see Eqn. (1161) . 

The next proposition shows that the maximal particle of the point process n* A 
has a shifted Gumbel distribution. 

Proposition 4. Assume that ip(X) = 0, A > 0, <p'(X) > 0. Let c(A) be as in 
Proposition fJJ Then, 

P[maxILj A <z]= exp | _ -y^ e ~ AZ | > z e M - 

2.3. Equilibrium states and their family structure. A point process 3 on R 
is called equilibrium state for the cluster distribution \ if the branching dynamics 
starting at 2 is stationary in time. With other words, if ~Y^t=\ is a realization of 
S and x , « G N, are copies of the point process \ ( au objects being independent), 
then the point process Y^tLi Tv t X nas the same law as S. Clearly, the point 
process IK is an equilibrium state for every A € K st , where K st = {AeR: <^(A) = 
0,Xip'(X) > 0}. Moreover, let 11* be the limiting point process for the branching 
dynamics started off the Poisson point process with intensity ce\, where c > and 
A G K st . Then, any mixture of the form 

(7) s(-)= / n* x (.)M<fc><*A)> 

where /i is a finite measure on the space (0,oo) x K st , is an equilibrium state. 
The next theorem describes the set of equilibrium states in the non-boundary case. 
A related result for systems without branching can be found [l^. Recall that a 
measure J on 1 is called non-lattice if it is not concentrated on an arithmetic 
progression of the form aZ + b, where a, b £ R. 

Theorem 2. Let x be a cluster distribution satisfying {T]) and having a non-lattice 
intensity measure J . Assume that the equation ip(X) = has two different solutions. 
Then, any equilibrium state 2 whose intensity measure is finite on bounded sets can 
be represented as a mixture of the form . 

Consider a branching dynamics {7r„ : n € No} starting off a point process JT* A . 

That is, ft n — SSi^V^Xn )i where ^vj is a realization of II£ and {x$ ■ 

n G No}, i G N, are copies of the branching random walk {y„ : n € No}, all objects 
being independent. Being stationary, the M-valued process {7r„ : n £ No} admits a 
two-sided stationary extension denoted by {7r„ : n € Z}. Let J 7 ^,! be the tr-algebra 
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generated by the random variables TTk{B), where B C M is a bounded interval and 
k = n, Ti—l, .... The cr-algebra of infinite past .F-oo is defined by J--00 = n„gz-?"<n. 

Proposition 5. Assume that ip(X) — 0, Ay'(A) > 0. Then, the a-algebra F-oo *s 
trivial, that is, every event in F-co has probability or 1. As a consequence, the 
process {TT n : n £ Z} is ergodic and mixing. 

Clearly, the union of two independent point processes distributed as n* ie and 
nx eA has the same law as LI* , ; , ci,C2 > 0. It follows that the point process 
IT* is infinitely divisible; see [18|, § 1 .61 . Moreover, the point process LI* is su- 



perposable in the following sense; see [17J ■ Recall that T u denotes a translation by 



it e R acting on the space M of discrete point configurations on R. Let IT and 
LT" denote two independent copies of II* . Then, for every u\,u<z,u € M such that 
e A«i _|_ gAu 2 _ e Aii we na ve the following equality in distribution: 

(8) T Ul Tl' + t U2 ti" =T u n* x . 

To see that |8]) holds, note that the Poisson point process with intensity e\ sat- 
isfies ([8]) and that superposability is preserved under clustering and taking weak 



limits. Superposable point processes have been characterized in 17] . As a conse- 
quence of [17] . we obtain 

Proposition 6. There is a unique in law point process T (depending on X) satis- 
fying maxL = a.s. and a constant c > such that ifTi, i € N, are independent 
copies of r and independently, * s a Poisson point process on R with in- 

tensity ce\, then 

00 

(9) n^^%,r, 

i=l 

It is convenient to think of © as of the decomposition of the point process IT* 



into the families T^T^, i 6 N; see 16J, § 1.6]. If we would continue the branching 
dynamics starting off IT* back into the past, then any two particles alive at time 
belong to the same family TyyA^i if and only if they have a common ancestor 
somewhere in the past. 

Proposition 7. Assume thatip(X) — 0, A > 0, (p'(X) > 0. If the cluster distribution 
X is critical, assume additionally that E[x(R) 2 ] is finite anrfP[x(R) = 1] < 1. Then, 
the point process T consists of finitely many particles a.s. if and only if the cluster 
distribution \ is subcritical. 

We have already noticed that even in the persistent branching dynamics the 
offsprings of any particle converge to — oo in a uniform way or die out. The next 
proposition says that a similar conclusion applies to the progeny of any family. 

(i) 

Let g n be the progeny of the family Ty/A^i at time n in the stationary branching 
dynamics starting off IT* . 

Proposition 8. Assume that <p(X) = 0, A > 0, <f'{\) > 0. If the cluster distri- 
bution x is critical, assume additionally that E[x(R) 2 ] is finite. If x * s critical or 
supercritical, then 

F[Vi E N : lim max^' = -00] = 1. 
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In particular, every family produces a.s. only finitely many leaders and finitely 
many particles above any fixed level. Also, it follows that there are a.s. infinitely 
many families (at time 0) which manage to produce a leading particle at some time 
in the sequel. Let us stress that the convergence in Proposition [8] is not uniform in 
i, since otherwise the whole branching system would drift to — oo like in the case 
of local extinction. 

Remark 2. If the cluster distribution \ allows only for integer- valued displacements 
of the particles, then it is natural to consider a branching dynamics on Z starting 
off a Poisson point process on Z with intensity e ~ Au< ^«- Our results have 

straightforward analogues in this settings. The non-lattice condition in Theorem [2] 
should be replaced by the assumption that J is aperiodic. 

2.4. The boundary case. We consider a system of particles starting off a Poisson 
point process on R with intensity e~ Xu du, where A > 0, and performing independent 
branching Brownian motions with drift — c, where c > 0. This means that each 
particle performs a standard Brownian motion with drift — c; after an exponential 
time it splits into two particles which behave in the same way as the original 
particle independently of each other. Branching Brownian motion is a continuous- 
time process, but its restriction to integer times belongs to the class of processes 
considered here. It is readily checked that (f(X) = (A 2 /2) — cA + 1. For c > \/2 the 
equation </?(A) = has two different solutions A± = c± Vc 2 — 2. By Theorem[T]thc 
measure e\ + corresponds to a persistent branching dynamics, while the measure 
e\_ corresponds to a branching dynamics which becomes locally extinct. In the 
boundary case c = y/2 there is only one solution A = y/2. Lalley and Sellke [HI 
conjectured that in this case the branching dynamics should be persistent. We will 
prove that this is not the case. 

Theorem 3. Consider a system of independent branching Brownian motions with 
drift — \/2 starting off the Poisson point process on R with intensity e Denote 
by 7r„ the point processes formed by the particles at time n. Then, ir n converges as 
n — > oo to the empty process. 

However, the second solution of = in the case c = V2 is not gone lost 
completely An easy computation shows that the signed measure v(du) — ue~^ u du 
is a solution of the convolution equation @ in the case c = It has been recently 
shown in [l[ that the branching dynamics starting off a Poisson point process with 
intensity —ue~^ u l u< odu converges to a non-trivial limiting point process. 

3. Proof of Theorem [T] 

3.1. Two lemmas on large deviations. Recall that I(z) = sup t6K (zi — (p(t)) is 
the Legendre-Fenchel transform of (p. We need the following simple Chernoff-type 
estimate. 

Lemma 1. Let {\n ■ n € No} be a branching random walk with cluster distribution 
X satisfying ([1]). Then, for every n £ N and a > mp'(0), 



(10) 



P[maxx„ > a] < e - nI{a/n) . 
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Proof. Let J™ = J * ... * J be the n-th convolution power of J. Note that is 
the intensity of \n- We have, for every t > 0, 

P[maxx„ > a] < ([a, oo)) < / e t{u - a) J [n] {du) = e -"^*-^)). 



The proof is completed by noting that sup t>0 (^t — <p(t)) — /(— ), where the supre- 
mum can be taken over t > since <p'(0). □ 

Lemma 2. Let {x n ■ n G No} & e a branching random walk with a supercritical 
cluster distribution \ satisfying (fTJ). If t > is such that <p'{t) < /3g, then 

lim -logE[e tmaxx "] = f3 t. 

n— too n 

Remark 3. If tp'(t) > /3o, then the limit is equal to ip(t); see Proposition [3] 

Proof. Write M n = maxx„. By (JSJ, we have limn^oo —M n = /3 a.s. on survival. 
It follows in a straightforward way that 

liminf-logE[e tMn l > f3 t. 

n~ VOO ft 

We prove the converse inequality. Trivially, E [e tMn lM n <p n 1 < e^ *". Fix e > 0. 
Then, by Lemma [l] 

E [e m "l M „>0 o n] < e* (z+£) "P[M„ > zn] < e etn ^ e (t*-/(«))n 

By convexity and since /(/So) = we have /(z) > I'{fio){z — /3q). Note also that 
i = < ^(^o) and ^ < 0. It follows that 

ze/3o+eN 

Bringing everything together, we obtain the statement of the lemma. □ 

3.2. Kallenberg's stability criterion. Our proof of the persistence part of The- 
orem[T]is based on the backward tree construction introduced by Kallenberg [l4| in 
the context of critical branching dynamics on R d starting off a homogeneous Pois- 
son point process. This approach has been generalized in § 2.4] (see also (20j ) 
to arbitrary branching dynamics on Polish spaces. 

Recall that we consider a branching dynamics on K d with spatially homogeneous 
cluster distribution \ having intensity J. We assume that A <E R d is such that 
<^(A) = 0. Let o be an individuum which is alive at time and is located at the 
origin in M. d . The idea of Kallenberg's method is to trace back the history of the 
individual o. First, we define the set of the direct ancestors of o. Let D be a 
probability measure on M. d defined by 

(11) D(B) = f e {X ' u) J{du), 

J B 

where B C K d is Borel. Since <p(\) = 0, D is indeed a probability measure. Let 
{£„ : n £ N}, be independent Revalued random variables with distribution D. Let 
S n = £i + • ■ • + £n, So — 0, be the corresponding random walk. It is convenient 
to think of — S n as of the position of the direct ancestor of o which is alive at time 
—n. So, — Si is the position of the father of o alive at time — 1, —52 is the position 
of the grandfather of o alive at time —2, etc. Now we define the set of relatives of o 



PERSISTENCE OF COMPETING SYSTEMS OF BRANCHING RANDOM WALKS 



9 



alive at time 0, i.e., the set of its brothers, cousins, etc. Let {x(z) : z £ ^. d } be the 
system of Palm measures of the point process Xi see [3, § 1.8] and [HI, Ch. 8,9]. 
Denote by {x\ z ) '■ z € the reduced Palm measures. Recall that x\ z ) ' IS obtained 
from X(z) by removing the point z. Conditioned on £i,£2,- •• let n n = Ylk=i 
rt G N, be independent point processes such that n n has the same law as X(£ )■ 
Think of —S n + z k , n , k = 1, . . . , k n , as of the locations of the brothers (alive at time 
— (n — 1)) of the ancestor of o located at — S n -i- Let also {xl*"'"' 1 '■ I € No}, where 
n £ N and k = 1, . . . , k n , be independent copies of the original branching random 
walk {x; : I £ N } generated by the cluster distribution x- Think of xP'™' as °f 
the point process describing the positions of the offsprings in the l-th generation of 
the individuum located at time — (n — 1) at — S n + Zk :n - Then, we define a point 
process p n by 

(12) ^=ET_ s „ +Zh: „(xi M 1 ) ). 

k=l 

Think of p n as of the point process describing the positions of the relatives of o 
which are alive at time and have the same ancestor at time — n as o. Let po = Sq. 



The next theorem is a specification of [16|, Proposition 2.4.3] to our setting 



Theorem 4 (Kallenberg's stability criterium). The branching dynamics starting 
off a Poisson point process with intensity e\ is persistent if and only if the measure 
p := X)^°=o P n * s l° ca My finite with probability 1. 

3.3. Proof of Part [T] of Theorem Q3 First we introduce some notation which is 
needed to handle the case d > 2 and may be ignored in the case d = 1. Without 
loss of generality we may assume that A = (1,0, ... ,0) for some I > 0, otherwise, 
we can rotate the coordinates. We denote by P : M. d — > R the projection operator 
P{1\, ■ ■ ■ , Id) — h- We may project the cluster distribution x to R via P. The 
resulting cluster distribution xp has intensity Jp{B) = J(P~ 1 B), where B C R 
is Borel, and the log-Laplace transform of Jp is given by <pp{t) — <f[t, 0, . . . ,0), 
t€R. Let Ip(z) = swp teR (zt — fp(t)) be the Legendre-Fenchel transform of (pp. 
We have <fp(l) — 0. Also, the assumption (A, Vy(A)) > implies that <p' P (l) > 0. 

By Thcorcm|4]it suffices to verify that p(A) is finite a.s. for A being a half-space 
set of the form [a, oo) x R rf_1 , a € R. Denote by A^ lK the a-algebra generated by 
the random variables {£„ : n £ N} and the point processes {k„ : n £ N}. We have, 
by (fT2"]) , Lemma [1] and the convexity of the function Ip , 



P[p n (A) > 0\A 6 , K ] 

k„ 

< £P [maxP(xi fe :" ) ) > P(S n ) - P(z k , n ) + a 

k=l 

/ , m fP(S n )-P(z k , n )+a 

i — i ^ 



*;=i 



V n - 1 



< g exp {-(„ - 1)Ip (£^>) + 7 P (|^) (P(z,„) - a)} 
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The expectation of £1 is given by E[£i] = J Rd uD(du) = V<p(X). By the law of large 
numbers, 

-> P(Vy>(A)) = ip'pil) > a.s. 

n 

By the properties of the Legendre-Fenchel transform, Ip(ip' p (l)) — lip' P (l) and 
I'p(<p'p(l)) — I- Take e > 0. It follows that with probability 1 we have for sufficiently 
large n, 

I p(^)>¥p(1)-£, l-s<I' P (^Pj<l + e. 
Thus, we have 

(13) nPn(A) > 0\At, K ] < Ce-^ l ^ n Y n , 
where {Y n : n € N} are random variables defined by 

fc =i jRd 

and A± e = (I ± e, 0, . . . , 0). We will show that E[F ra ] is finite. Define a function 
/ : M — >• [0, +oo] (where M is the space of all locally finite counting measures on 
R d ) by f(f3) = j Rd {e {X+ " u) J r& {X -" u) )P{du). Then, 

E[Y n }= [ E[f( X { u) )}J(du)<[ E[f( X(u) )}J(du)=E[f(x)}-nx(^% 

where the last step follows from the definition of the Palm measure. It follows that 

(14) E[Y n ] = (e 1 ^^ + e^ A -^)e^°). 

Recall that in Part[T]of Theorem[T]we assume that lf' P {l) > 0. It follows from (fT5|) 
and ((HI) that 

OO OO 

J2^[pn(A) > 0] = ^E[P[p„(A) > 0\Az, K }} < oo. 

n— 1 n— 1 

By the Borel-Cantelli lemma, p(A) — Y^=o Pn{A) ls finite a.s. By Theorem 2] the 
proof is completed. 

3.4. Proof of Part [2] of Theorem [TJ 

Proof of Proposition^ Recall the definition of 7r n in ([6]). The point process 7r* :— 
Si^i V+maxx (i) * s a P°i sson point process on M. with intensity E[e Amaxx "] -e~ Xu du. 
For every f3 > (3$ we have 

OO OO OO 

(15) ^P[max7r n > j3n] < ^ E[<((£n, oo))] = -J2(e' Xfln E[e Xmaxx -}) < oo, 

n— 1 n— 1 n— 1 

where the last step follows from Lemma [5] (note that <p'(X) < (3o). By the Borel- 
Cantelli lemma, we have limsup^^^ ^ max7r„ < /3q. The converse inequality fol- 
lows immediately from ([5]). □ 
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Proof of Part[^ of Theorem^ Let first d = 1. For concreteness, assume that A > 

0. Then, <p'(0) < 0, /3o < and Proposition [2] implies that lim I1 _ i . 00 max n n = — oo 
a.s. It follows that tt„ converges weakly to the empty process. 

We consider the case d > 2. Without restriction of generality let A = (1,0, ... ,0) 
with I > 0. Let A C R d be a bounded Borel set. We will show that with probability 

1, tt„(A) = for sufficiently large n. Let E n =Rx [-n 2 , n 2 ] d - 1 and E n = R d \E n . 
Let N n (respectively, N n ) be a random variable counting the number of particles 
of the point process 7r„ which are in the set A and whose ancestor at time was 
in the set E n (respectively, in E n ). Clearly, it n (A) — N n + N n . After elementary 
transformations we obtain 



E[N n ]= / e-( x ^jW(A-u)du = / e-^D^(a- E n )da, 

JE n J A 



where = D * ... * D is the n-th convolution power of the probability measure D 
defined in (fTTj) . By the Cramer-ChernofT bound there is c > such that D^(a — 
E n ) = 0{e- cn ) uniformly in a £ A. Hence, V[N n > 0] < E[N n ] = 0{e~ cn ). By the 
Borel-Cantelli lemma, N n = for all but finitely many n a.s. 

Let us show that with probability 1, N n — for all but finitely many n. Consider 
the progeny of all particles which start at time in the cylinder set R x [0, 
Let us project these particles onto the first coordinate R via the projection operator 
P; see Section [5751 We obtain a branching dynamics {iri^ : n G No} on R starting 
off a Poisson point process with intensity e~ lx . By Proposition [2] (see also (jX5j) ) 
there are P < 0, c> such that P[max7ri P) > j3n] = 0{e~ cn ). Divide the set E n 
into (2n 2 ) d ~ 1 sets which are translates of the set R x [0, and apply to each of 
them the above considerations. Let R n be the largest among the first coordinates of 
the offsprings in the n-th generation of the particles which started in E n . It follows 
that F[R n > /3n] = 0(n 2d ' 2 e- cn ) = 0(eT cn l 2 ). By the Borel-Cantelli lemma, 
R n — > — oo a.s. It follows that with probability 1, N n = for all but finitely many 
n. This completes the proof. □ 



3.5. Proof of Theorem [3j Let Xn be the point process formed by the particles at 
time n in a single branching Brownian motion with drift —y2. Write M n = max% n . 
It suffices to show that 



(16) lim E[e 



n— f oo 



Fix e > 0. Note that the expected number of particles at time n in the branching 
Brownian motion is e n . Denoting by $ the tail of the standard normal distribution 
function we have a trivial estimate 



P [M n >y]< E[ Xn ([y, oo))] = e"<i f V + ^ | < C 
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It follows that for a sufficiently large c > 0, 



E 



,V2M r , 



M n >[cy/n\ 



<c eV ^ kp t M « ^ fc ] 

fc=[c % /n] 

c 



< c 

< e. 



e 2 (is 



Write l n — ^= logn. By Proposition 3 of Bramson [5[ there is B = B(c) such that 
for all — l n < y < c*/n, we have the following estimate 

B 



' [M n >y]< 



n 3/2 e V2y 



This is stated in 5] for c = 1 only, but the proof is valid for every c > 0. It follows 
that 



E 



e V2M„ i 

e J— [i n ]<M„<[cVn] 



[cy/n\ 

< e^ {k+1) P[M n >k} = 

fc=-[i»] 



Trivially, we have lim„_ ; . 0O E 
we obtain flTBl) . 



- J-Af n <-!„ 



= 0. Bringing everything together, 



4. Proof of Propositions [H El El [S 

4.1. Proof of Propositions [3] and [4j Assume first that tp(t) = and > /?o> 
which, in fact, implies that t(p'(t) > 0. Then, the branching dynamics {7r„ : n £ No} 
starting off a Poisson point process with intensity &t is persistent by Theorem [T] 
In particular, it follows that lim„_ i . 00 P[7T n ((z, +oo)) = 0] exists in (0, 1) for at least 
one z £ R. On the other hand, 

(17) P[max7r„ <z] = F[ir n ((z, +oo)) = 0] = exp j-i E [ e * max *"] ■ e~ t2 

Hence, c(i) = lim^oo E[e tmaxx ™] exists in (0, oo). 

If ip(t) and <p'(t) > /3o, we consider another branching random walk {x„ : 
n £ No} obtained by adding to {\ n ■ n £ No} a drift —(p(t)n/t, that is x„ = 
T- v (t) n /t{Xn)- The log-Laplace transform of xi satisfies (p{t) = and <j?'(i) > (3q — 
= Po, where /3q is the analogue of /3q for the function I(z) = sup sgR (sz — <p(s)). 
By the above, 

lim e -¥>(*)»E[e* roaxx ™] = lim E[e* max *"] G (0,oo). 

This completes the proof of the first part of Proposition [3] The second part follows 
from Lemma [5J Proposition @] follows by taking the limit n — > oo in (fTT)) . 
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4.2. Proof of Proposition [H The idea of the proof is as follows. At time 0, 
the density of particles at —un is e Xun , and the probability that a particle at —un 
will produce at time n an offspring, say, on the positive half-axis is approximately 
e -nj(u). cf. Lemma [T] Thus, the ancestor at time of the particle which is the 
leader at time n is most likely to be found at position approximately —un, where 
u is such that Xu — I(u) is maximal. The maximizer is given by u = <f'(X), whence 
the statement. 

Let us provide a rigorous argument justifying this. Fix some e > and write 
a±(e) = <p'(X) ± £. We estimate the probability -P^"(e) that some offspring of some 
particle which started at time below the level — a + (e)n will be above the level 
—y/n at time n: 



where we have used a substitution u = nv + yfn and Lemma [TJ The last inequality 
holds if n is sufficiently large. The function g(v) := Xv — I(v) is concave and we have 
g(<p'(X)) — 0, g'(ip'(X)) = 0. The global maximum of g is attained at v = <f'(X). 
By concavity there are c\ > 0, c 2 > such that g(v) < — c\ — c 2 (w — a + (e/2)) for 
all v > a + (e/2). Hence, 



Analogous considerations apply to P~(s), the probability that some offspring of 
some particle which started at time above the level — a_(e)n will be above the 
level — yfn at time n. Hence, Y^=iPn{ e ) < 00 an d Y^=i Pn{ £ ) < 00 • By the 
Borel-Cantelli lemma, only finitely often a particle starting not in the interval 
(— a+(e)n, — a_(e)n) will be able to produce a particle above the level —\ph at 
time n. Also, we have 



where the last inequality holds since E[e AmaxXl1 ] remains bounded by Proposition^ 
Hence, only finitely often the leading particle will be below the level —-Jn. The 
statement of the proposition follows. 

4.3. Proof of Propositions [7| and [HI We will use the notation introduced in 
Section l3~2l If the cluster distribution is critical, we always assume that E[x(K) 2 ] 
is finite. 

Lemma 3. If the cluster distribution \ is subcritical or x(^) = 1 a - s -> then p(M) < 
oo a.s. Otherwise, p(R) = oo a.s. 




P 



max (U l + maxXn ) > -\fn 



_Ui<-u + (e)n 
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Proof. In the subcritical case we have E[%(R)] < 1 and hence, 

oo oo oo 

E[p(R)} = ]TEMR)] = 1 + ^E[fc n ]E[x„-ipR)] = 1 + ^^(EIxCK)])"- 1 < oo. 

ri—0 n— 1 n— 1 

It follows that p(M) < oo a.s. In the non-subcritical case we write 

P[p n (K) > 0] > ¥[k n > l.x^W > 0] = P[fc„ > 1] • P[x„-i(R) > 0]. 

If we assume that P[x(R) = 1] < 1, then P[fc„ > 1] > is a constant independent 
of n. Further, in the critical case we can find c > such that for all n G N, 
P[x„-i(R) > 0] > c/n; see p. 19]. (Here the finiteness of E[x(R) 2 ] is needed). 
In the supercritical case, we even have P[x»-i(R) > 0] > c. By the Borel-Cantelli 
lemma, p(M) = oo a.s. □ 



Lemma 4. Let A C M be a Borel set such that if ft £ A, then T u /3 <E A for every 
u G M. Let 93(A) = 0, A > 0, <£>'(A) > 0. TTien, P[T G A] = i/ and onfy i/ 
P[p G A] = 

Proof. As already observed in Section I2.3[ the point process IT* is infinitely divis- 
ible. Denote by \I/ its canonical (Levy) measure; see [18|, Ch. 2]. Note that ^> is 
a generally infinite measure on M. On the one hand, the cluster representation of 
IT* given in Proposition [6] implies that we have 



(A) = c J P[T„r G A]e- Xu du. 



Thus, P[r e A] = if and only if = 0. On the other hand, by Theorem 2.4.4 

of [16| we have 

9 (u) {A)=F[T u peA], ueR, 

where {^(V) : u G R} is the family of Palm measures of *f>. Thus, P[p e 4] = 
if and only if \I/( U )(A) = for all u 6 R. Finally, by the definition of the Palm 
measure we have for every Borel set Bel, 



J B 

Thus, * (u) (A) = for all u G M if and only if =0. □ 

Proof of Proposition [?J The proof follows from Lemma [3] and Lemma |4] by taking 
A to be the set of all point configurations f3 G M such that /3(R) = 00 (in the 
subcritical case) or /3(R) < 00 (otherwise). □ 

Lemma 5. Let (p(X) — 0, A > 0, <^'(A) > 0. If x is critical or supercritical, then 
J R e ru T(du) < 00 for every r > rg, where ro > is the smallest solution of the 
equation (p(r) — ip'(\)r. 

Proof. By Lemma H we need to show that J R e ru p(du) < 00 a.s. for all r > r$. 
Recall that is the a-algebra generated by {£„ : n G N}. Denote by Jr u \ the 
intensity of the point process X(«)> w e M. Fix e > 0. We have S n /n — > iff (A) by 
the law of large numbers. It follows that with probability 1, S n /n > f'(X) — e for 
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sufficiently large n. Hence, 

(18) E [ f e ru p{du)\A^\ = 1 + f] [ e ru T_ Sn {J [n ~ 1] * %„))(<M 

Jm J n=1 Jr 

OO p 

n=l jR 

OO 

n=l 

oo 

< B + e v ^ 52 e(~ r(p '^ + ' p ( r ) +re ^ n Y n , 

where Y n — J R e ru J^ n ^(du) and B is an a.s. finite ^-measurable random variable. 
Let us show that E[Y n ] is finite. Recall that the distribution of is e Xz J(dz). 
Using the properties of the Palm measures, we have 



Wn] =1(1 e ru J iz) (du)) e Xz J(dz) 

JR \JR / 

= [e( f e ru+Xz X(z)(du)) J(dz) 



= E / / e ru+Xz x(du) X (dz) 
Jr Jr 

— e <p(r)+<p(\) 

So, M[Y n ] is finite. If e > is chosen such that e < (fi'(X) — -<p(r), then the series 
on the right-hand side of (IT51) converges a.s. It follows that E[L e ru p{du)\A{\ < oo 
a.s. Hence, J R e ru p(du) < oo a.s. □ 

We are now in position to prove Proposition |SJ Let 5p i be a realization of 

the point process T from Proposition [51 Consider a branching dynamics {7„ : n 6 
No} starting off T. That is, let {\n^ ■ n € No}, i £ I, be independent copies of the 
branching random walk {y„ : n G No}. Then, j n — ^2 ieI Tp i (xn^)- We need to 
show that 

(19) lim max7„ = — oo a.s. 

n— >oo 

Let _4r be the a-algebra generated by T. We have 

F[Pi +maxx { n ] > -en\A r ] < e^" 5 "^ < e -™'(-<V'(-*)A. 

We have J(— e) > and I'{—e) > ro if e > is sufficiently small since lim^o /(— e) = 
7(0) = inf tS R<^(i) < and lim e _>.o e ) = -f'(O) = argmint^ > ro- It follows from 
Lemma [5] that 

oo 

P [max7„ > — ert|^r] < oo a.s. 

n=l 

An application of the Borel-Cantelli lemma completes the proof of (QjJJ . 
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5. Proof of Theorem [2] and Proposition [5] 

The proof is based on the general results of [la ]. By Theorem 2.3.6 of [l6| a 
measure v on M. d which is finite on bounded sets is an intensity of some equilibrium 
state if and only if the branching dynamics starting off a Poisson point process with 
intensity v converges weakly to a point process with intensity v (such measures v 
are called stable). Thus, by Theorem [TJ the intensity v of an equilibrium state S is 
a linear combination with non- negative coefficients of the measures of the form e\, 
where A € K st . 



By Theorem 4.7.28 in 16[ any equilibrium state S with intensity v is a mixture 
of the form (JT)) if the following two conditions (called aperiodicity and regularity) 
are satisfied: 

(1) For every probability measure a on K which is absolutely continuous with 
respect to the Lebesgue measure we have 

(20) lim \\a * - a * L>[" +1 ]|| = 0, 

n— >oo 

where D is the probability measure defined in (fTTj) and || • || denotes the 
total variation. 

(2) For every bounded Borel set A C K and every e > 0, 

(21) lim v{z e K : J [n] (A - z) > e} = 0. 

n— >oo 

See Proposition 3.8.4 in [l6| for the first condition and Theorem 4.7.19 in [l6| 
for the second one. The aperiodicity condition holds as long as the measure J 
is non-lattice by Theorem 11.10.4 and Proposition 11.9.1 of 18]. Let us verify 
the regularity condition. In the subcritical case the measure on the left-hand side 
of (|2"Tj) is for sufficiently large n. In the non-subcritical case, we need a different 
argument. Without restriction of generality let v — e\, where ip(X) = 0, A > 0, 
<p'(\) > 0. It suffices to consider the case A — (a, oo), where net, We have, for 
every t > 0, 

J [n] (A- z) < e - t[a - z) [ e tu J [n] (du) = e -*(«-*)+"*(*). 

JR 

It follows that 

(22) {zeR: J [n] {A-z) > e} C + ^ + a,oo 

We can choose t > such that ip(t) < 0. For example, take t = A — 6 for sufficiently 
small S > 0. The e^-measure of the set on the right-hand side of (|22p goes to as 
n — > oo. This completes the proof of Theorem [2] Proposition [5] follows from the 
regularity condition (|2"Tj) by Theorem 4.7.19 of [161 ] . 
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